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Abstract. The reduction of dynamical systems has a rich history, with many important apphcations related 
to stability, control and verification. Reduction of nonlinear systems is typically performed in an "exact" 
manner — as is the case with mechanical systems with symmetry — which, unfortunately, limits the type of 
systems to which it can be applied. The goal of this paper is to consider a more general form of reduction, 
termed approximate reduction, in order to extend the class of systems that can be reduced. Using notions 
related to incremental stability, we give conditions on when a dynamical system can be projected to a lower 
dimensional space while providing hard bounds on the induced errors, i.e., when it is behaviorally similar to 
a dynamical system on a lower dimensional space. These concepts are illustrated on a series of examples. 



1. Introduction 

Modeling is an essential part of many engineering disciplines and often a key ingredient for successful designs. 
Although it is widely recognized that models are only approximate descriptions of reality, their value lies 
precisely on the ability to describe, within certain bounds, the modeled phenomena. In this paper we consider 
modeling of closed-loop nonlinear control systems, i.e., differential equations, with the purpose of simplifying 
the analysis of these systems. The goal of this paper is to reduce the dimensionality of the differential 
equations being analyzed while providing hard bounds on the introduced errors. One promising application 
of these techniques is to the verification of hybrid systems, which is currently constrained by the complexity 
of high dimensional differential equations. 

Reducing differential equations — and in particular mechanical systems — is a subject with a long and rich 
history. The first form of reduction was discovered by Routh in the 1860's; over the years, geometric reduction 
has become an academic field in itself. One begins with a differential equation with certain symmetries, i.e., 
it is invariant under the action of a Lie group on the phase space. Using these symmetries, one can reduce 
the dimensionality of the phase space (by "dividing" out by the symmetry group) and define a corresponding 
differential equation on this reduced phase space. The main result of geometric reduction is that one can 
understand the behavior of the full-order system in terms of the behavior of the reduced system and vice 
versa |MW74I IvdSSlI IBKMM96] . While this form of "exact" reduction is very elegant, the class of systems 
for which this procedure can be applied is actually quite small. This indicates the need for a form of reduction 
that is applicable to a wider class of systems and, while not being exact, is "close enough". 

In systems theory, reduced order modeling has also been extensively studied under the name of model re- 
duction |BDG96| lASGOO] . Contrary to model reduction where approximation is measured using norms 
we are interested in L°° norms. The guarantees provided by L°° norms are more natural when applica- 
tions to safety verification are of interest. More recent work considered the exact reduction of control sys- 
tems |vdS04l ITP04] based on the notion of bisimulation which was later generalized to approximate bisimu- 
lation IGP05I ITabOTI IGPOT] . 

We develop our results in the framework of incremental stability and our main result is in the spirit of existing 
stability results for cascade systems that proliferate the Input-to-State Stability (ISS) literature. See for 
example [Son06' and the references therein. A preliminary version of our results appeared in the conference 
paper [TAJP06j '. 
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2. Preliminaries 

A continuous function 7 : Rq , is said to belong to class /Coo if it is strictly increasing, 7(0) — and 

7(r) — > 00 as r — > 00. A continuous function (3 : M.^ x is said to belong to class /C£ if, for each fixed 

s, the map /3(r, s) belongs to class /Coo with respect to r and, for each fixed r, the map /3(r, s) is decreasing 
with respect to s and P{r,s) ^0 as s ^ 00. 

For a smooth function : R" R™ we denote by Tel) the tangent map to </> and by T^j^ the tangent map to 
<j) at X G R". We will say that is a submersion at a: € R" if Tx4> is surjective and that T(j) is a submersion if 
it is a submersion at every x £ R". When is a submersion we will also use the notation ker(T0) to denote 
the distribution: 

ker(T0) = {X : R" ^ R" \ Tc/) ■ X ^ 0}. 
The Lie bracket of vector fields X and Y will be denotes by [X, Y] . 

Given a point x E R", |a;| will denote the usual Euclidean norm while ||/|| will denote esssup^gjo t] 1/(^)1 foi' 
any given function / : [0, t] ^ R", t E R+. 

2.1. Dynamical and control systems. In this paper we shall restrict our attention to dynamical and control 
systems defined on Euclidean spaces. 

Definition 2.1. A vector field is a pair (R",X) where X is a smooth map X : R" R". A smooth curve 
x( • , : I ^ R", defined on an open subset / of R including the origin, is said to be a trajectory of (R", X) 
if the following two conditions hold: 

(1) x(0,x) ^x; 

(2) ^x(t,x) = X{x{t,x)) for all t E I. 

A control system can be seen as an under-determined vector field. 

Definition 2.2. A control system is a triple (R",R™,i^) where F is a smooth map F : R" x R™ ^ R". A 
smooth curve Xu(-, x) : I ^ R", defined on an open subset J of R including the origin, is said to be a trajectory 
of (R", R™, F) if there exists a smooth curve u : / ^ R™ such that the following two conditions hold: 

(1) Xu(0,a;) = x; 

(2) ^Xu(t,x) = F(Xu(i,a;),u(t)) for all t E I. 

We have defined trajectories based on smooth input curves mainly for simplicity since the presented results 
hold under weaker regularity assumptions. 

3. Exact reduction 

For some dynamical systems described by a vector field X on R" it is possible to replace X by a vector 
field Y describing the dynamics of the system on a lower dimensional space, R™, while retaining much of the 
information in X. When this is the case we say that X can be reduced to Y. This idea of (exact) reduction 
is captured by the notion of (/)-related vector fields. 

Definition 3.1. Let : R" R™ be a smooth map. The vector field (R", X) is said to be (fy-related to the 
vector field (R™, Y) if for every x E R" we have: 

(3.1) T,4>-X{x)^Y o4>{x). 



The following proposition, proved in jAMR88j . characterizes (/)-related vector fields in terms of their trajecto- 
ries. 
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Proposition 3.2. The vector field (M",X) is (j>-related to the vector field {R"\Y) for some smooth map 
(j) : M" M™ iff for every x G M" we have: 

(3.2) 0ox(t,x) =y(tXa;)), 

where ii(t,x) and y(t,y) are the trajectories of vector fields X and Y, respectively. 

For (/)-related vector fields, we can replace the study of trajectories x(-,a;) with the study of trajectories 
y(-,(/)(x)) living on the lower dimension space E™. In particular, formal verification of X can be performed 
on Y whenever the relevant sets describing the verification problem can also be reduced to M™. 

If a vector field and a submersion (j) are given we can use the following result, proved in |MSVS85] . to determine 
the existence of c/i-related vector fields. 

Proposition 3.3. Let be a vector field and let <j) '■ E" — s- M™ he a smooth submersion. There exists a 

vector field {W"\Y) that is (j)-related to iff: 

[ker{T(t)),X] C ker{T(j)). 



In an attempt to enlarge the class of vector fields that can be reduced we introduce, in the next section, an 
approximate notion of reduction. 



4. Approximate Reduction 



The generalization of Definition 3.1 proposed in this section requires a decomposition of K" of the form 
M" = M™ X M*^ . Associated with this decomposition are the canonical projections tt™ : M" and tt/c : 

M" R'' taking 3 x ^ {y, z) e M™ x to 7r™(a;) = y and nkix) = z, respectively. 

Definition 4.1. The vector field (M",X) is said to be approximately TT„i-related to the vector field (M™,y) 
if there exist a class /Coo function 7 and a constant c G Rq such that the following estimate holds for every 
X G K": 

(4.1) |7r„ o x(t, x) - y(t, 7r„(x))| < 7(|7rfc(a;)|) + c. 



Note that when X and Y are TTm-related we have: 

|7r„ ox(i,a;) - y(i, 7r„(a;))| = 0, 



which implies (4.11. Definition 4.1 can thus be seen as a generalization of exact reduction captured by 
Definition |3.1[ Similar ideas have been used in the context of approximate notions of equivalence for control 
systems GP07 . 

Although the bound on the gap between the projection of the original trajectory x and the trajectory y of 
the approximate reduced system is a function of a;, in concrete applications the initial conditions are typically 
restricted to a bounded set of interest. The following result has interesting implications in these situations. 

Proposition 4.2. If{R^,X) is approximately TTm-related to {R"\Y) then for any compact set S C M" there 
exists a S G M"*" such that for all x G S the following estimate holds: 

(4.2) |7r™ox(t,a;) -y(t,7r™(a;))| < (5. 



Proof: Let 5 = max^jgc 7(|7rfe(x)|) + c. The scalar 6 is well defined since 7(|7rfc( • )|) + c is a continuous map 
and C is compact. □ 

From a practical point of view, approximate reduction is only a useful concept if it admits characterizations 
that are simple to check. In order to derive such characterizations we need to review several notions of 
incremental stability. 
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4.1. Incremental stability. In this subsection we review two notions of incremental stability which will be 
fundamental in proving the main contribution of this paper. We follow |BMOO] and |Ang02| . 

Definition 4.3. A control system (M", K™, F) is said to be incrementally uniformly bounded-input-bounded- 
state stable (lUBIBSS) if there exist two class /Coo functions 71 and 72 such that for each xi,X2 € M" and for 
each pair of smooth curves Ui, U2 : / ^ M™ the following estimate holds: 

(4.3) \x^^{t,xi) - yi^^{t,X2)\ < 7i(|xi - a:;2|) + 72(||ui - U2II) 

for all t e /. 



In general it is difficult to establish lUBIBSS directly. A sufficient condition is given by the existence of an 
IUBIBSS Lyapunov function. Note, however, that lUBIBSS only implies the existence of a lUBIBSS Lyapunov 
function with very weak regularity conditions [BMOOJ. 

Definition 4.4. A function V : E" x M" M.^ is said to be an lUBIBSS Lyapunov function for control 
system (M", M™, F) if there exist a ^ S and class /Coo functions a, a, and n such that for every xi,X2 G K" 
and ui,U2 G M™ the following holds: 

(1) Ixi - X2I > C =^ a{\xi - X2\) <V{xi,x2) <a{\xi ~ X2\); 

(2) fj.{r) > r + ^ for r e M+; 

(3) IX1-X2I >m(|wi-M2|) ^ V<0. 



A stronger notion than lUBIBSS is incremental input-to-state stability. 

Definition 4.5. A control system (M^jR"*, F) is said to be incrementally input-to-state stable (IISS) if there 
exist a class JCC function /? and a class /Coo function 7 such that for each xi,X2 € K" and for each pair of 
smooth curves Ui,U2 : I ^ M™ the following estimate holds: 

(4.4) \x^^{t,xi) - yi^^{t,X2)\ < f3{\xi - X2\,t)+j{\\ui - U2II) 



Since /3 is a decreasing function of t we immediately see that (4.4 1 implies (4.3 1 with 71 (r) = /3(r, 0) and 
72(?') — i{t), r e MJ. Once again, IISS is implied by the existence of an IISS Lyapunov function. See |Ang02| 
for a converse result when the inputs take values in a compact set. 

Definition 4.6. A function F : K" x K" ^ K+ is said to be an IISS Lyapunov function for the control 
system (M",M™,F) if there exist class /Coo functions a, a, a, and ji such that: 

(1) a{\xi -X2\) < V{xi,X2) < a{\xi ~ X2I); 

(2) |xi - X2I > At(|Mi - U2I) =^ V" < -a(|a;i - X2I). 



4.2. Fiberwise stability. In addition to incremental stability we will also need a notion of partial practical 
stability. 

Definition 4.7. A vector field (M",X) is said to be fiberwise practically stable with respect to TTfc if there 
exist a class /Coo function 7 and a constant c G M.q such that the following estimate holds: 

||7rfc(x(-,a;))|| < 7(|7rfe(a;) |) -|- c. 
Fiberwise practical stability can be checked with the help of the following result. 

Lemma 4.8. A vector field (R.",X) is fiberwise practically stable with respect to tt^ if there exist two /Coo 
functions, a and a, a constant d € M,^, and a function V : M" M such that for every x E M" satisfying 
|7rfe(a;)| > d we have: 



(1) a{\Trk{x)\)<V{x)<a{\TTk{x)\), 

(2) V<0. 
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4.3. Existence of approximate reductions. In this subsection we prove the main resuh providing sufficient 
conditions for the existence of approximate reductions. 

Theorem 4.9. Let (W^,X) be a fiberwise practically stable vector field with respect to tt^ and let F — TiTm-X : 
X M'^ ^ M", viewed as a control system with state space W"\ be lUBIBSS. Then, the vector field (M",r) 
defined by: 

for every y G M'" is approximately Tr,n-related to (M",X). 

Proof: By assumption, control system (M'",E'^. F = Ttt^ o X) is lUBIBSS. If we denote by y a trajectory of 
F we have: 

|yvi(t,j/i) - yv2(^,y2)| < ii{\yi - ^2!) + 72(l|vi - V2II). 

In particular, we can take: 

yi ^ y2 ^ T^mix), Vi = TTfe O x(-,x), V2 = 0, 

to get: 

|7r„ o x(i,a::) - y(i,7r„(a;))| = \y7TkOx{t,x){t,TTm{x)) - yo{t,TTjn{x))\ 

= |yvi(i,7r™(a;)) -yo(i, 7r„(a;))| 
< 72(||vi||) = 72(||7rfe ox(-,a;)||). 

But it follows from fiber practical stability of X with respect to tt^ that: 

IjTTfe o x(-,a;)|| < 7(|7ri.(a;)|) + c. 

We thus have: 

|7r„ox(t, x) -y(t, 7r™(x))| < 72(7(|7rfe(a;)|) + c) 

< 72(Ai7(|7rfc(a;)|)) +72(A2c), 

for some constants Ai,A2 S M.^ . This concludes the proof since 72(Ai7(| • |)) is a class /Coo function and 
72(A2c) G K+. □ 

Theorem |4.9| shows that sufficient conditions for approximate reduction can be given in terms of ISS-like 
Lyapunov functions and how reduced system can be constructed. Before illustrating Theorem 4.9 with several 
examples in the next section we present an important corollary. 



Corollary 4.10. Let and be vector fields satisfying the assumptions of Theorem 4-9 Then, 

for any compact set S C M" there exists a S > such that for any x Cz S and y G iTmiS) the following estimate 
holds: 

\-Kjn o x(t,x) - y{t,y)\ < 5 

Proof: Using the same proof as for Theorem |4.9[ except picking yi — nm{x) and y2 = y, it follows that: 

|7r,„ ox(i,a::) - y{t,y)\ < 7i(|7r„(x) - 2/|) + 72(Ai7(|7rA;(a;)|)) +72(A2c). 
The bound S is now given by: 

."^^^ 7i(|7r,^(x) -y|) + 72(Ai7(|7rfe(a::)|)) +72(A2c) 
which is well defined since S x TTm{S) is compact. □ 





Figure 2. A trajectory of the full order system (red) vs. a trajectory for the reduced system 
(blue) for R — 5, 10, 20, 40 (from left to right and top to bottom, respectively). 



5. Examples 

In this section, we consider examples that illustrate the usefulness of approximate reduction. 
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Example 5.1. As a first example we consider the ball in a rotating hoop with friction, as described in Chapter 
2 of |MR99] and displayed in Figure [l] For this example, there are the following parameters: 

m = mass of the ball, 

R = Radius of the hoop, 

g = acceleration due to gravity, 

/i — friction constant for the ball. 



uj - 

m 



sin 9 cos — ^ sin 



The equations of motion are given by: 

LU 

(5.1) 9 
where 9 is the angular position of the ball and lu is its angular velocity. 

If TT^^ : ^ M is the projection ^^^{uj, 9) = ui, then according to Proposition 3.3 there exists no vector field Y 
on M which is TTt^-related to X (as defined by (5.1 1). However, we will show that Y{uj) — T(^^ q-^tt^ ■ X{uj, 0) is 
approximate Tr^j-related to X. 

First, we use: 



V 



1 



mR?LO^ + mgR{l — cos ( 



1 



mi?^^2 sin^ ( 



as a Lyapunov function to show that (5.1 ) is stable. Note that V{lj, 9) = for {uj, 9) = (0, 0) and V{uj, 9) > Q 
for {uj,9) 7^ (0,0) provided that R^^ < g, which we assume. Computing the time derivative of V we obtain: 

V = -^iR'^Lo^ < 0, 



thus showing stability of (5.1 1. We now consider a compact set C invariant under the dynamics and restrict 
our analysis to initial conditions in this set. Such a set can be constructed, for example, by taking {x G 
I V{x) < c} for some positive constant c. Note that stability of (5.1 1 implies fiberwise stability on C since 
TTjniC) is compact. 



To apply Theorem 4.9 we only need to show that: 



9 



X{uj, 9) = -—u) + sm9cos9 - —sm( 
m R 



is lUBIBSS on C with 9 seen as an input. We will conclude lUBIBSS by proving the stronger property of 
IISS. Consider the function: 

U = -{UJI - UJ2T- 



Its time derivative is given by: 

U 



(Wl - UJ2) 



{uji — UJ2) + sin 6*1 cos 9i 



R 



sin 9i — £^ sin ^2 cos 6*2 + — sin ^2 



< 



< 



- — (^1-^2) + \ui-ijJ2\ 

m 



(5.2) 



R 

m 
2m 



sin 9i — sin 6*2 cos 6*2 



R 

sin 9i cos 9i 
. 9 



R 



smt/2 



{uJi - ^2)^ + \iLJi - UJ2\L\(- 



{UJI - LU2y 



JJ_ 

2m 



{loi 



Iwi - Ld2\L\9l 



where the second inequality follows from the fact that sin 9 cos 9 — sin 6* is a smooth function defined on 
the compact set ng(C) and is thus globally Lipschitz on ng(C) (since its derivative is continuous and thus 
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R 



Figure 3. A graphical representation of the pendulum on a cart mounted to a spring. 



bounded on any compact convex set containing TTg(C)) with Lipschitz constant L. We now note that the 
condition: 

|wi - W2I > \0l - 02] 

A* 



makes the second term in (5.2) negative from which we conclude the following implication: 



\UJI -UJ2\ > \0l -6*2! 



Zm 



showing that U is an IISS Lyapunov function. We can thus reduce (5.1 1 to: 



-Ul. 



m 



Projected trajectories of the full-order system as compared with trajectories of the reduced system can be 
seen in Figure [2] here ^ = m = 1 and ^ = 0.1. Note that as i? ^ 00, the reduced system converges to the 
full-order system (or the full-order system effectively becomes decoupled). 



Example 5.2. We now consider a pendulum attached to a cart which is mounted to a spring (see Figure |3|. 
For this example, there are the following parameters: 



M = 


mass of the cart, 


m = 


mass of the pendulum. 


R = 


length of the rod, 


k = 


spring stiffness. 


9 = 


acceleration due to gravity, 


d = 


friction constant for the cart. 


b = 


friction constant for the pendulum 
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Figure 4. A projected trajectory of the full-order system (red) and a trajectory for the 
reduced system (blue) for d = 0.001, 0.01, 0.1, 1 (from left to right and top to bottom, respec- 
tively). 



The equations of motion are given by: 



V 



M + m 



^ — (mRuj^ sin 6 + mq sin 9 cos 6 — kx — dv + — cos 9\ 

nsm^0\ R J 

— mRu)^ sin cos 



(5.3) 



i?(M-f TO sin^ 
— {m + M)g sin + kx cos + dv cos 















TO 





where x is the position of the cart, v its velocity, is the angular position of the pendulum and uj its angular 
velocity. 



If TT 



is the projection 7r(^x_y-^{x,0,v,uj) — {x,v) and X is the vector field as defined in (5.3 1, the 



goal is to reduce X to R by eliminating the and uj variables and thus obtaining Y defined by: 



Y{x, v) = Ti^x,e,v,uj)T^{x,v) ■ X{x, 0, 0) 



jjidv + kx) 
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The objective is now to show that X and Y are approximately Tr^jj^-j-related. In particular, note that the 
reduced system, Y , is linear while the full-order system, X, is very nonlinear. This will be discussed in more 
detail after proving that they are in fact approximately related. 

Stability of X, and in particular fiberwise stability, can be proven as in the previous example by noting that 
X is Hamiltonian for d = = and using the Hamiltonian as a Lyapunov function V . Consider now the 
control system: 

(5.4) F{{x,v),{e,Lo)) = T7r^,,,yX{x,0,v,Lo) 



1 (r 

n sin^ 6 V 



M + m sin^ 9 



niRuj sin 9 — k 



X 



-mg sin 9 cos 9 — dv -\ cos ( 

R 



with 9 and w regarded as inputs. To show that F is lUBIBSS we first rewrite (5.4 1 in the form: 

F iix, v), (9, u))) — (mRu)^ sin 9 ~ kx ~ dv ~ mRuj cos 9^ 

and consider the following IISS candidate Lyapunov function: 



Its time derivative is given by: 



mR 



(^Lul sin 9i — ui cos 9i — sin 92 + 1^2 cos (I'l — ^2)- 



m + M 

Using an argument similar to the one used for the previous example, we conclude that: 

'2jjiRF 

\vi - V2\ > I (6*1,^1,101) - (6*2,^2,^^2)1, 

with L the Lipschitz constant of the function cj^ sin 9 — uj cos 9, implies: 

^^-2(;dni/)(^^-^^)'' 

thus showing that X is IISS and in particular also lUBIBSS. That is, we have established that X and Y are 
approximately 7r(a:.u)-related. 

In order to illustrate some of the interesting implications of approximate reduction, we will compare the 
reduced system, Y, and the full-order system, X, in the case when R = m = k = b= l and M — 2. It follows 
that the equations of motion for the reduced system are given by the linear system: 

X \ ( \ \ ( X 



so we can completely characterize the dynamics of the reduced system: every solution spirals into the origin. 
This is in stark contrast to the dynamics of X (see (5.3 1) which are very complex. The fact that X and Y are 
approximately related, and more specifically Theorem |4.9[ allows us to understand the dynamics of X through 
the simple dynamics of Y. To be more specific, because the distance between the projected trajectories of X 
and the trajectories of Y is bounded, we know that the projected trajectories of X will "essentially" be spirals. 
Moreover, the friction constant d will directly affect the rate of convergence of these spirals. Examples of this 
can be seen in Figure [4] where d is varied to affect the convergence of the reduced system, and hence the full 
order system. 
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